Cooling rate, heating rate and aging effects in glassy water 
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We report a molecular dynamics simulation study of the properties of the potential energy land- 
scape sampled by a system of water molecules during the process of generating a glass by cooling, 
and during the process of regenerating the equilibrium liquid by heating the glass. We study the 
dependence of these processes on the cooling/heating rates as well as on the role of aging (the time 
elapsed in the glass state). We compare the properties of the potential energy landscape sampled 
during these processes with the corresponding properties sampled in the liquid equilibrium state 
to elucidate under which conditions glass configurations can be associated with equilibrium liquid 
configurations. 

PACS numbers: 



One recent activity in the physics of supercooled liq- 
uids and glassesPjQ is the search for the conditions un- 
der which a glass can be considered a liquid whose struc- 
tural properties have been "frozen" during the prepa- 
ration process. If the glass can be connected to a liq- 
uid state, then a thermodynamic description of the glass 
state can be developed|3,0,IE0- Many routes can bring 
a system to an arrested disordered state ItJ, such as va- 
por deposition, pressure induced amorphization, hyper- 
quenching or standard cooling. Only the last two provide 
a continuous path from the liquid to the glass state and 
hence are the best candidates for studying the connection 
between glass and liquid configurations. 

Since water can be glassified by cooling only using hy- 
perquenching techniques (i.e. with rates of the order of 
10^ K/sQ), understanding the connection between the 
liquid state and glasses generated with different cooling 
rates is important. When the hyperquenched glass of wa- 
ter is properly annealed at T = 130 kQ, a reproducible 
weak endothermic transition is observed which has been 
associated with the calorimetric glass transition temper- 
ature Tg. For water, the standard cooling rate is not suf- 
ficiently fast to overcome crystallization, so active debate 
|9l ll0lll]| concerns how to relate the Tg of hyperquenched 
water to the unmeasurable Tg of the slowly cooled glass. 

In this work we use molecular dynamics (MD) simu- 
lations to address the relation between liquid and glass 
configurations. We reproduce the same procedure fol- 
lowed experimentally to generate glasses. We use two 
cooling rates qc differing by almost three orders of mag- 
nitude. We also study the heating rate and aging effects 
in the glass, an aspect novel to numerical simulations. 
We work in the framework of the potential energy land- 
scape (PEL) approach, in which the GA^-dimensional con- 
figurational space — defined by the 2>N center of mass 
coordinates and by the 3iV Euler angles — is partitioned 
into a set of basins, each associated with a different local 



minimum of the PEL [13. Il3l |. We focus in particular 
on the depth of the local minimum closest to the system 
point, and on the local curvatures of the PEL around 
the local minimum. We discover that with slow cooling 
rates, the glass retains a configuration very similar to 
a configuration sampled by the liquid at higher T, and 
hence all structural properties of the glass can be related 
to the structural properties of the liquid. In this case, 
a Active temperature can be defined, the temperature at 
which the glass configuration is sampled by the liquid. 
In the case of a fast cooling rate, aging phenomena are 
very active already during the cooling process. Dynamics 
moves the configuration in regions of the PEL which are 
never explored in equilibrium p^ . Further aging at low 
temperature increases the differences between the glass 
and the liquid. When this is the case, the glass does 
not possess a structure corresponding to the equilibrium 
liquid at any temperature and hence it is not possible to 
associate an unique Active T with the glass configuration. 

We perform MD simulations for a system of iV = 216 
molecules at fixed density, p = 1 g/cm'^, interacting with 
the SPC/E potentialll5|, with periodic boundary condi- 
tions. Interactions are cut off at a distance of r = 2.5cr 
(cr parameterizes the Lennard- Jones part of the SPC/E 
potential) and reaction field corrections are added to ac- 
count for the long range interactions. Quantities are aver- 
aged over 32 independent trajectories. We perform three 
types of MD calculations: (i) cooling scans at constant 
rate, starting from equilibrium liquid configurations at 
T = 300K, (ii) heating scans at constant rate (from glass 
configurations at « 5 K) , and (iii) aging runs at constant 
aging temperature (at Tagc = 100 K and Tagc = 180 K, 
where significant aging effects are observed). The two 
cooling rates are Qc — —3 x 10^° K/s and qc ~ —10^^ K/s, 
and the two heating rates are qh — +3 x 10^° K/s and 
qh — -I-IO^'^ K/s. We will denote by fast-quenched glass 
the glass obtained with the fast cooling rate and by 
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slow-quench glass the glass obtained with the slow cool- 
ing rate. An averaged slow scan requires a simulation 
lasting 320 ns, close to the maximum possible by our 
method. This limitation prevent us from studying larger 
systems . The location of the system on the PEL is 
studied performing numerical minimizations (conjugate 
gradient algorithm) along the runs to estimate the clos- 
est local minimum configuration or inherent structure 
(IS), its energy e/s and the set of 6N eigenvalues ujf 
of the Hessian matrix^^. The local curvatures around 
the minimum define in the harmonic approximation the 
multi-dimensional parabolic shape of each PEL basin. As 
a global indicator of shape we use the shape function de- 
fined as 

^'^ = N ^> 

where Aq = I kJ/mol. At density p = Ig/cm^ the tem- 
perature dependence of ejs a,nd Sjs for the SPC/E model 
has been previously calculated |19l|. 

Focusing on landscape properties has the advantage 
that vibrational features are suppressed, and hence the 
remaining temperature dependence reflects structural 
changes. The landscape approach is particularly useful 
for comparing properties of different configurations. In- 
deed, at constant volume, in equilibrium, all landscape 
properties (i.e., e/s, Sjs) are functions only of T. By 
eliminating T, we can establish a direct relation between 
Sis and ejs, i.e. between properties of the landscape 
sampled in equilibrium. A set of configurations with the 
same e/s can be considered representative of equilibrium 
configurations if and only if Sjs satisfies the equilibrium 
Sisi^is) landscape relation. We compare the relation 
between shape and depth during cooling scans, heating 
scans and aging with the equilibrium relation. This al- 
lows us to test if the region of configuration space ex- 
plored during the scans is identical to the region explored 
by the equilibrium liquid at some higher temperature. 

Figure^ shows eis{T) for the fast and the slow cool- 
ing rates, together with equilibrium values from Ref . [l^ . 
The equilibrium data show that in the liquid state, the 
system explores deeper and deeper basins on lowering T 
and that at each T, basins with a well-defined energy 
depth are populated. The fact that different depths are 
associated with different T provides a mechanism for as- 
sociating a configuration (whose e/s value is known) to 
a fictive temperature Tf , as suggested in Ref. 0| . In the 
fast cooling case the liquid loses equilibrium at a very 
high temperature (around 290K), i.e., when the system 
is still sampling basins of high energy. In the slowest cool- 
ing case, the system loses equilibrium around 220 K. For 
the fast cooling case, the range of T in which e/s changes 
is significantly larger, and more energy (from the T where 
equilibrium is lost to the final T) is released during the 
cooling scan. 

Figure^ shows e/s for the fast and slow heating rates 
starting both from the fast-quenched glass and the slow- 
quenched glass. For the fast heating rate, to reach equi- 



librium T must reach a value larger than the T at which 
the same system lost equilibrium on cooling (Fig. ^i). 
In the slow heating case, a new phenomenon appears for 
the fast-quenched glass, the system starts to explore PEL 
basins with lower and lower e/s as it relaxes toward equi- 
librium, significantly changing the location of the system 
in configuration space. 

Figures andCJi show the 5/s(e/s) relation for the 
scans reported in Figs.^ and^D, respectively. One sees 
that in the slow cooling case, the cooling scan proceeds 
via a sequence of equilibrium configurations (since the 
'5/s(e/s) relation coincides with the equilibrium one), i.e. 
basins with same depth explored during the slow scan and 
those in equilibrium have also the same shape. However, 
fast cooling generates a glass (the black square) which we 
see lies outside the equilibrium curve of Fig. QJ;. In this 
case, basins explored during fast cooling that have the 
same depth as those found in equilibrium have smaller 
curvatures. During the slow heating process <S/s(e/s) 
moves even further from the equilibrium line (Fig. ^) . 

We now focus on the trajectory of the fast-quenched 
glass during aging at two different aging temperatures 
Tage = lOOif and Tagc = 180K. For both rates we find 
that e/s is a monotonically decreasing function of aging 
time iago (Fig. |2i), suggesting that the system changes 
basin by finding specific directions on the PEL which in- 
volve low energy barriers. For long iago, e/s ~ log(iagc), 
as found in atomic liquids 113, Il3l- During aging, two dif- 
ferent behaviors are again observed for the PEL <S/s(e/s) 
relation (Fig.|2|D). In the high Tage scan, the system ap- 
proaches equilibrium configurations, while for low Tago 
the opposite is seen. Note (Fig. ^p) that aging at low 
Tage moves the system even further from the equilibrium 
5/s(e/s) curve. 

Next we show (Fig.|3l heating scans for different aging 
times after isothermal aging at Tago = 100 K. Data for 
Tage — 180 K (not shown) are qualitatively similar. The 
IS energy of the starting configuration is lower the longer 
the time spent at Tage = 100 K before beginning to heat. 
We see (Fig. |3|;) that e/s is monotonic on heating. Fig- 
ure|3l3 shows the results for the fast heating rate. In this 
case, T increases so fast that the system does not have 
time to find deeper IS. Therefore, there are no aging ef- 
fects and e/s grows monotonically with T. From Figs.|2t 
andlJb, we see that for the fast heating rate, equilibrium 
is reached at a much higher temperature compared to 
the case of the slow heating rate, in full agreement with 
the no-aging (tage = 0) observed in Fig. The 

aging effects on 5/s(e/s) are shown in Figs.|3; (for slow) 
and |3i (for fast) heating rates, confirming that in the 
slow heating case the fast-quenched glass keeps explor- 
ing new basins never seen in equilibrium before meeting 
the equilibrium curve. 

In summary, we have shown that the fast-quenched 
glass significantly differs from the slow-quenched glass. 
To the extent that this model captures key features of 
real water, our work provides support to the idea that, 
if crystallization would not interrupt the cooling process 
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in water, the slow-quenched glass would be character- 
ized by a different Tg compared to the fast-quenched 
glass 0, 0|- The analysis of simulated configura- 
tions has allowed us to clarify the differences between 
the slow- and the fast-quenched glass. The differences 
lie not only in the expected difference in e/5 but, more 
significantly, in the fact that the fast-quenched glass dur- 
ing cooling loses contact with the liquid state and starts 
to explore regions of the PEL which are never explored 
in equilibrium. The fast-quenched glass does not pos- 
sess a structure corresponding to the equilibrium liquid 
structure at any temperature, and hence it is not pos- 



sible to associate with it a unique Active temperature. 
Finally, we have found that aging in the glass phase pro- 
vides a mechanism for reducing differences between the 
liquid and the glass structures only if the temperature at 
which the system ages is sufficiently high. Remarkably, 
at low Tago, differences between the liquid and the glass 
structures appear to increase with tage. 
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FIG. 1: (a) Cooling of equilibrium liquid configurations at 
T = 300 K to T « K with a fast {q^ = -lO" K/s) and slow 
{qc = —3 X 10^'' K/s) cooling rate to generate fast-quenched 
and slow-quenched glasses, respectively. While during fast 
cooling the system is out of equilibrium already at T = 290 K, 
slow cooling allows the system to sample the equilibrium IS 




FIG. 2: (a) Time dependence of eis and (b) relation between 
basin shape Sis and basin depth during aging of the fast- 
quenched glass at two temperatures Tago — 100 K and Tage = 
180 K. For both cases, for long tage, eis ~ log(tagc) (dashed 
lines). Note that during aging at low Tagc, the system explores 
basins for which the relation between shape and depth does 
not correspond to the equilibrium Sis{eis) relation. At large 
Tage, aging moves the system closer to equilibrium. At each 
Tage, open diamonds correspond to (from right to left) tage = 
20 pS, tage = 200 pS, tage = 2 ns, and tage = 20 ns. 
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FIG. 3: Inherent structure energy eis and basin shape func- 
tion Sis during heating of the fast-quenched glass with (a,c) 
slow and (b,d) fast heating rates after aging for different times 



tagc at a temperature Taj 



100 K. Results for T=180 K are 



qualitatively similar. Filled circles correspond to data for equi- 
librium liquid configurations from Rcf. [19]For comparison, in 
(a) and (b) we show eis during the heating scan from the 
slow-quenched glass. In (c) and (d), we only show Sis in the 
heating scans after aging for tage = 20 ps and tage = 20 ns. 



